Vortex in a trapped Bose-Einstein condensate with dipole-dipole interactions 
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We calculate the critical rotation frequency at which a vortex state becomes energetically favorable 
over the vortex-free ground state in a harmonically trapped Bose-Einstein condensate whose atoms 
have dipole-dipole interactions as well as the usual s-wave contact interactions. In the Thomas- 
Fermi (hydro dynamic) regime, dipolar condensates in oblate cylindrical traps (with the dipoles 
aligned along the axis of symmetry of the trap) tend to have lower critical rotation frequencies than 
their purely s-wave contact interaction counterparts. The converse is true for dipolar condensates 
in prolate traps. Quadrupole excitations and centre of mass motion are also briefly discussed as 
possible competing mechanisms to a vortex as means by which superfluids with partially attractive 
interactions might carry angular momentum. 
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I. INTRODUCTION 

The achievement of Bose-Einstein condensation (BEC) 
in a trapped gas of 52 Cr atoms by the Stuttgart group 
PJ is the first instance of a condensate with large dipole- 
dipole interactions. In comparison to alkali atoms, which 
have a maximum magnetic dipole moment of fj, = 
(fiB is the Bohr magneton), chromium atoms possess 
an anomalously large magnetic dipole moment of fj, = 
6^b- The long-range part of the interaction between two 
dipoles separated by r, and aligned by an external field 
along a unit vector e (in this paper we shall always as- 
sume the dipoles are aligned by an external field), is given 
by 



^dd(r) = - — e t ej — - — 
47T r 



(1) 



where the coupling Cdd = MoM 2 depends on the square of 
the dipole moment. Chromium atoms also have shorter 
range isotropic interactions which are asymptotically of 
the van der Waals-type and so fall off as r~ 6 . At ultra- 
low temperatures the de Broglie wavelength of the atoms 
is much larger than the range of the isotropic interactions 
which can consequently be handled within perturbation 
theory by using the usual delta-function pseudo-potential 
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U(r) = 4irh 2 a s 5(r)/m = gd(r) 



(2) 



characterized solely by the s-wave scattering length a s . 
A measure of the strength of the long-range dipole-dipole 
interaction relative to the s-wave scattering energy is 
given by the dimensionless quantity 



£dd 



Cdd 

~3g" 



(3) 



For 87 Rb one finds e dd « 0.007, and Na has e dd ~ 0.004. 
Taking the s-wave scattering length of 52 Cr to be a s = 



105ae [1|, one obtains £ dd ~ 0.144. In a further signif- 
icant development, the Stuttgart group has also studied 
Feshbach scattering resonances between 52 Cr atoms 

3]- 

This means that the magnitude and sign of a s , and thus 
also the value of e dd , can be controlled. 

The principal effect of the anisotropy of the dipole- 
dipole interactions upon a stationary, trapped, conden- 
sate will be to distort its aspect ratio so that it is elon- 
gated along the direction of the external field 0, @ . This 
is in contrast to the case of purely isotropic interactions, 
for which the aspect ratio of the condensate matches that 
of the trap. By adopting an elongated aspect ratio along 
the direction of polarization, the condensate achieves a 
lower energy by placing more dipoles end-to-end in which 
configuration they are attractive, and reducing the num- 
ber of repulsive side-by-side interactions. When e dd ex- 
ceeds a certain value, which in general depends on the as- 
pect ratio the trap as well as the total number of atoms 
N, mean- field theory predicts that the condensate be- 
comes unstable to collapse 0, H, El El • The partially 
attractive nature of the dipole-dipole interaction is also 
res pon sible for introducing a curious 'roton' minimum 
13, llOj l into the Bogoliubov excitation spectrum of a uni- 
form dipolar BEC, reminiscent of that found in the ex- 
citation spectrum of liquid helium II. It is possible that 
this roton minimum is indicative of an instability towards 
the formation of a density wave 0, 0] . 

The fact that the interaction between dipoles aligned 
by an external field along the axis of symmetry of a cylin- 
drical trap are on average repulsive in an oblate (pancake 
shaped) BEC, but, conversely, are on average attractive 
in a prolate (cigar shaped) BEC, means that the sign of 
the dipolar mean-field energy can be controlled via the 
aspect ratio of the trap. In this paper we are interested 
in calculating the effect that this change of aspect ra- 
tio of the trap has upon the vortex state in a trapped 
dipolar BEC. The case of a vortex in a trapped BEC 
with repulsive short-range interactions, as described by 
Eq. with a s > 0, has been extensively studied both 
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theoretically, e.g. 0, 0], and experimentally, e.g. 
[S [H Eg. Attractive short -range interactions have, 
on the other hand, received comparatively less attention, 
although a rotating BEC with attractive interactions is 
of considerable fundamental interest E3- ln particular, 
when the interactions are attractive it is intuitively plau- 
sible that vortex formation will be suppressed because it 
is energetically expensive to form the vortex core since 
this requires moving atoms from the centre of the BEC 
where they interact with many other atoms and placing 
them on the edges where the interact with fewer [20| . 
Instead, it has been proposed that the angular momen- 
tum might be preferentially absorbed into other types of 
excitation, such as centre of mass motion or shape os- 
cillations, which cause much less disturbance to the con- 
densate density profile and internal correlations |19l |2l| . 
We also note that a recent theoretical study of a two- 
dimensional BEC with purely attractive short-range in- 
teractions found a new class of vortices in the form of 
bright ring solitons ■ When it comes to dipolar inter- 
actions, which are partly attractive, partly repulsive, one 
might therefore anticipate rotational properties similar 
to those of standard BECs with either attractive or re- 
pulsive interactions, depending upon whether the dipolar 
BEC is prolate or oblate, respectively. 

The latest experiments on chromium j23| have detected 
the first evidence of dipolar interactions upon the dynam- 
ics of a BEC by observing the ballistic expansion of the 
BEC when the trap is switched off. The results are in 
good agreement with theoretical predictions |24j giving 
us some confidence in the theory of dipolar BECs which 
has been built up over the last six years. In the following 
we will discuss a vortex in a rotating dipolar BEC in the 
Thomas-Fermi (hydrodynamic) regime, which is defined 
as being when the quantum zero-point kinetic energy due 
to the confinement by the trap is negligible in compar- 
ison to the trapping and interaction energies. This is 
the relevant regime when there are a large number of 
atoms. We shall make extensive use of the exact results 
for dipolar BECs in the Thomas-Fermi limit reported in 
references [2^ and [2(j . The key insight of those papers 
was that the dipolar mean-field potential that the atomic 
dipoles feel, and which is non-local, i.e. depends on the 
entire atom distribution, can be expressed in terms of 
derivatives of a scalar potential that is a solution of the 
Poisson equation with the dipolar density as a source. In 
other words, this scalar potential is formally equivalent 
to the electrostatic potential one would get if the den- 
sity distribution were not one of dipoles but of charges, 
and thus standard mathematical techniques known from 
electrostatics can be applied. Three recent theory papers 
[27Ll2§ . l2^| have dealt with the structure of vortex lattices 
in dipolar BECs, but our aim here is to investigate the 
differences between dipolar and non-dipolar condensates 
in the case of a single vortex. Our main concern will 
be the calculation of the critical rotation frequency f2 c 
necessary to make a single vortex energetically favorable 
over the ground state. Note that reference [2!j also deals 



with the single vortex case, but their focus was somewhat 
different. One of their principle results was that in the 
case of dipolar interactions plus attractive contact inter- 
actions (negative scattering length) the vortices develop 
a 'craterlike' structure which is probably connected with 
the density wave instability mentioned above. We shall 
come back to this point in Section ^ however we shall 
not explicitly consider negative scattering lengths here 
because we wish to use some exact results that only hold 
in the Thomas- Fermi limit. A dipolar BEC with attrac- 
tive short-range interactions is unstable to collapse in the 
Thomas- Fermi limit and so we limit ourselves to a s > 0. 



II. ENERGY FUNCTIONAL FOR A DIPOLAR 
BEC WITH A VORTEX 

Consider a condensate at T — with a wave func- 
tion ^(r) normalized to the total number of atoms 
J \^f(r)\ 2 d 3 r = N. The general form of \&(r) when there 
is one singly quantized vortex present can be expressed in 
cylindrical coordinates (p, (j>, z) as 4'(r) = \^(p, z) \ expi0 
|2j. This state carries an angular momentum L z = N%. 
In a frame rotating at angular velocity £1 about the z- 
axis the energy of the condensate becomes E' = E — VtL z , 
where the angular momentum L z and energy E are those 
pertaining to the laboratory frame. However, creating a 
vortex costs energy. Denoting the energy of the BEC in 
its ground state without any vortex (L z = 0) by Eg and 
the extra energy needed to make a single vortex by E v , 
we can write the energy of the vortex state in the rotating 
frame as 

E'=E +E V - flL z . (4) 

Thus we see that for a BEC at equilibrium in a trap ro- 
tating at angular velocity £1 the vortex state becomes en- 
ergetically favorable when f2 exceeds a critical rotational 
velocity 



To obtain f2 c we thus need to evaluate the extra energy 
associated with the formation of the vortex. The total 
energy functional for a trapped dipolar BEC can be writ- 
ten as 

E to t = -Ekinotic + -Strap + E sw + -^dd (6) 

where, in terms of the condensate wave function ^(r), 

Akinetic = J d 3 r **(r)V 2 *(r) (7) 

is the kinetic energy, and 

Strap = ^|d 3 r|vI/(r)| 2 [p 2 +7 V] (8) 
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is the energy due to the harmonic trap Vtrap 
(m/2)w 2 [p 2 + 7 2 ^ 2 ], where p 2 = x 2 



■ y 2 and 



7 = lu z /u x 



(9) 



is the ratio of the trap frequencies. Note that in this pa- 
per we shall assume that both the trap and the conden- 
sate are cylindrically symmetric about the z direction, so 
in particular lo x = lu v , and the external field responsible 
for aligning the dipoles is along the z axis. 

The total mean-field interaction energy can to a good 
approximation [j| be written as the sum of two parts. 
The first is due to the isotropic short-range interactions 
which give rise to pure s-wave scattering 



E sw = | J d 3 r|*(r) 



(10) 



and the second is due to dipole-dipole interactions 

E dd=\[ d 3 r dV|*(r)| 2 C/ dd (r - r')|*(r')| 2 . (11) 

The long-range and anisotropic nature of the dipole- 
dipole interaction, Udd( r ), makes the calculation of E^d 
a non-trivial exercise. However, in the Thomas-Fermi 
regime some exact results (within mean-field theory) are 
available. In particular, an exact solution of the vortex- 
free Thomas-Fermi problem for a dipolar BEC in a har- 
monic tra p y ields a density profile which is an inverted 
parabola [23, H|| . This solution is similar to the famil- 
iar non-dipolar case 0, except that the aspect ratio of 
the cloud in the dipolar case is no longer identical to 
that of the trap but is stretched along the direction of 
the polarizing field, as mentioned above. By applying 
scaling transformations to this exact solution one can 
also describe the lowest energy collective excitations of 
the condensate [2q|. as well as the expansion dynamics 
when the trap is turned off. In the Thomas- Fermi regime 
the critical ratio of the dipolar and s-wave interactions 
above which instabilities occur is £dd = 1 30] . The oc- 
currence of these instabilities is indicated by the appear- 
ance of imaginary frequencies for the collective excita- 
tions. Another particularly novel instability in a trapped 
dipolar BEC, which originates from the long-range na- 
ture of the interatomic interactions, is the appearance of 
a local minimum in the mean-field potential outside the 
condensate when e d d > 1- This stimulates the formation 
of a 'Saturn-ring' around the condensate as atoms tunnel 
into the minimum p^ |. 

As explained at length in references [25| and j2||, the 
dipole-dipole interaction can be re-expressed so as to take 
advantage of the large technical machinery that exists to 
deal with electrostatic interactions. In general one can 
write 



E ^^\j d3 r »(r)*dd(r) 



where $dd( r ) is the non-local dipolar mean- field potential 



$ dd (r) = J dV U dd (r - r')n(r') 



(13) 



For dipoles aligned along z this dipolar mean-field po- 
tential can in turn be expressed in terms of a fictious 
electrostatic potential 0(r) 



$dd(r) 



-C, 



<id 



d 2 . n(r) 
3^ (r) + 3 



(14) 



which is obtained from the "charge" distribution n(r) in 
the usual way 



, n(v') 



r - r' 



(15) 



where, of course, n(r) and (j>(r) satisfy Poisson's equation 
V 2 = — n(r). This formulation of the problem allows us 
to immediately identify some generic features of dipolar 
gases. For example, if n(r) and hence </>(r) are uniform 
along the polarization direction (z) then the non-local 
part of the dipolar interaction vanishes because of the 
d 2 <j>/dz 2 operator. A very prolate pencil-like condensate, 
i.e. one which has R z 3> R x ,R y will therefore be largely 
unaffected by the non-local part of the interaction except 
perhaps at the very ends of the condensate if the den- 
sity profile has a large curvature there. Thus, in the very 
prolate (or entirely homogenous) case the contribution of 
dipolar interactions to the mean-field potential reduces to 
that of the local term <f>dd( r ) — > — Cdd^(r)/3 whose effect 
is then simply to modify the magnitude of the local mean- 
field potential arising from the pure s-wave interactions. 
A related example to a very prolate BEC is a BEC with 
a generic smooth outer density profile but which is, how- 
ever, penetrated by a pencil-like vortex or vortices: the 
dipolar potential would be expected to peak around the 
ends of each vortex where the density profile is rapidly 
changing in the z-direction. One might conjecture that 
the curious 'craterlike' density profile around the ends of 
vortices in dipolar BECs that has recently been seen in 
simulations [23 is connected to this effect. 

Let us consider the case of a parabolic density profile 
of the form 



n bg (r) = n ( 1 



P_ 

Rl 



z 

m 



(16) 



As stated above, a parabolic density profile is an exact so- 
lution of the hydrodynamic equations for a non-rotating 
BEC in the Thomas-Fermi limit even in the presence of 
dipolar interactions. Here the subscript "bg" stands for 
"background" , i.e. for a condensate without a vortex but 
on top of which we will superimpose a vortex later on. 
One finds that the dipolar mean-field potential inside a 
condensate with the parabolic density profile ljl6|l is given 
by EHH 



(12) ***(p, z ) 



bg/ n-oCdd 



Rl 



2z 2 ( 3 p 2 -2z 2 



In this expression 
/(«) 



2 + k 2 [4-3S(k)] 
2(1 — K 2 ) 



(18) 



(17) 
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is a function of k which is the aspect ratio of the BEC, 



Rx/Rz 



(19) 



and monotonically decreases from its value of / = 1 at 
k = 0, passing through zero at k = 1, tending towards 
/ = —2 as k — > oo. The function S(«) upon which /(k) 
depends is itself dependent upon whether the condensate 
density profile is prolate, in which case 



1 



\/l — K 2 1 — y/l — K 
or oblate 



, i + VT - ^ , -. / i + n 

In . lor k < 1 (prolate) 



(20) 



arctan vkM for k > 1 (oblate) 



(21) 

The value of the aspect ratio k is give n by the solution 
of a transcendental equation 0, 124. l25l l2q 



o 2 

3k e d d 



1 



- 1 



+ (e dd - 1) (k 2 - 7 2 



= 0. 



(22) 

Note that according to Eq. 1|17[) . the dipolar mean- field 
potential $^(p,z) is in general saddle-shaped [2^1 , and 
thus inherits the anisotropic partially attractive/partially 
repulsive character of the dipolar interactions. The sad- 
dle shaped mean-field potential causes an elongation of 
the BEC along the polarization direction which can be 
viewed as type of magnetostriction. This is illustrated 
for five different trap aspect ratios in Fig. ^ The total 
energy (E tlap + -^sw + E^d) associated with the vortex- 
free Thomas-Fermi solution l|16|l can be calculated to be 
M 



En 



N 
14 



muj 2 R 2 
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15 N 2 g 
28tt RIR Z 



[1 - £dd/(«)] 



(23) 

The quantity Eq is the ground state energy referred in 
Eq. Q. 

Having now summarized what is already known con- 
cerning the mean-field potential inside a non-rotating 
dipolar BEC, let us turn to the form of the density pro- 
file of a rotating BEC with a vortex. We shall adopt the 
following variational ansatz for the density profile of an 
iV-atom condensate with a single vortex 



n(r) 



|*(r)|< 
n I 1 



£1 _ 1 
Rl R 2 

n bg (r) +n v (r) 



1 



/3 2 



(24) 
(25) 



where we have observed that the ansatz can be written 
as the sum of two terms: the background Thomas-Fermi 
parabolic profile Hb g (r) as already given in Equation 11 tip 
and the vortex profile n v (r) 



n v (r) 



-n 



2 



p 2 +p 2 



Rl R 2 



(26) 
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FIG. 1: Aspect ratio k — R x /R z of a vortex-free dipolar 
BEC as a function of scattering length a s (measured in units 
of the Bohr radius ao) in the Thomas- Fermi limit. In de- 
scending order, the curves are for traps with aspect ratios 
7 = Lo z juj x = 6, 3, 1, 0.4, 0.2. When a s — > oo the s-wave con- 
tact interactions dominate and the aspect ratio matches that 
of the trap, i.e. k = 7, given by the right hand asymptote 
of each curve. When a s — > the dipolar interactions dom- 
inate and their magnetostrictive effect reduces k below the 
pure s-wave value. In the limit a s — > a dipolar BEC will 
collapse towards a chain of end-to-end dipoles and mean-field 
theory will breakdown. In this figure Cdd = Mo(6mb) 2 and the 
minimum value of scattering length shown is o 3 = 17.5ao. 



Ansatz (|24f) assumes that the vortex, whose core size is 
parameterized by 0, is superimposed on a background 
parabolic density profile with radii R x — R y and R z . 
Note that the variational ansatz l|24ll has the correct p 2e 
dependence of the density as p — > (with £ = 1), which it 
must have in order to satisfy the Gross-Pitaevskii equa- 
tion for a vortex of I circulation quanta [2l| . It has also 
got the correct asymptotic form (in terms of even powers 
of l/p) that the solution must have for p — > 00. The 
central density no is fixed by normalization to be 



"0 



8irR 2 x R Zl 



1 + ^-p 2 + 5/5 4 



5p 2 (1 + P 2 f' 2 arctanh[l/Jl + (3 2 \ ) (27) 



where 



P = &/R* 



(28) 



is the ratio of the vortex core size to the transverse radius 
of the BEC. 

We proceed by minimizing E to t as a function of the 
three variational parameters R x , R z and (3. The terms 
^kinetic, -^trap and E sw can all be calculated analytically 
in a straightforward albeit laborious manner. The results 
are presented in Appendix A. To evaluate the dipole- 
dipole energy functional we begin from Equation l|ll|) 
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and substitute in the density ansatz 125(1 



E dd - - j d 3 r dV n bg (r)f/ dd (r - r> bg (r') 
+ i J d 3 r dV n bg (r)£/ dd (r - r')n v (r') 
+ i / d 3 r dV n v (r)(7 dd (r - r> bg (r') 
+ i y d 3 r dV n v (r)L/ dd (r - r> v (r'). 

Noting that the two cross terms between the vortex and 
background densities are identical, i.e. the integral is in- 
variant under exchange of the coordinates r and r', the 
dipolar energy functional can be written 

E dd = i J d 3 r dV n hg (r)U dd (r - r> bg (r') 
+ J d 3 r d 3 r' n v (r)U dd (r - r')n bg (r') 
+ i J d 3 r dV n v (r)L/ dd (r-r> v (r'). 

We shall now restrict ourselves to situations where the 
size of the vortex core is much smaller than the radius 
of the condensate, i.e. j3 -C 1. This is consistent with 
the spirit of the Thomas-Fermi approximation for the 
background parabolic envelope, and should hold for con- 
densates containing a large number of atoms, providing 
the system does not become very prolate. To this end we 
shall drop the vortex-vortex part of the dipolar energy 
functional since this has an extra factor of p 2 in compar- 
ison to the cross term. The remaining two terms can be 
expressed as 

E dd « \j d 3 r dVn bg (r)C/ dd (r-r')n bg (r') 
+ j d 3 r dV n v {v)U dd (v - r')n bg (r') 
= \J d 3 rn bg (r)$^(r) + J d 3 r n v (r)^(r!|29) 

The first term is just the dipolar energy of the conden- 
sate without a vortex, whilst the second gives the dipo- 
lar 'interaction' between the vortex and the background 
density profile. Both terms have been reduced to sin- 
gle integrals of their respective density profiles over the 
known quadratic function $ dd (/?, z), as given by Equa- 
tion l)17p. which is the mean- field potential due to the 
background inverted parabola density profile. These in- 
tegrals can therefore be evaluated relatively easily and 
the results are given in the Appendix. 



III. RESULTS OF MINIMIZATION OF THE 
ENERGY FUNCTIONAL 

A. Oblate trap 

To illustrate the oblate case we shall consider 150,000 
52 Cr atoms in a trap having frequencies oj x — u> y — 
2tt x 200 rad/s, lj z — 2n x 1000 rad/s, so that 7 = 5. 
The harmonic oscillator length of the trap along the x- 
direction is then a bo = v^7 muJ x = 0.986/im. We numer- 
ically minimize the total energy functional -Etot with re- 
spect to {/3, R x , k} using the Mathematica routine NMin- 
imize. In order to bring out the dependence of the vari- 
ous quantities upon the relative strength of the two types 
of interactions (s-wave and dipolar) we plot the results 
as a function of the s-wave scattering length a B , since 
this quantity can be controlled in an experiment [3|. As 
a s — > the dipolar interactions dominate, whereas when 
a s — > oo the dipolar interactions become insignificant. 
The magnitude of the magnetic dipole interaction is here 
assumed to be fixed at the appropriate value for 52 Cr, 
namely C dd = /io(6/is) 2 . Note that in the case of dipoles 
induced by electric fields one has C dd = E 2 a 2 /eo |3lj . 
where a is the polarizability and E is the electric field 
strength, so the size of the electric dipolar interaction 
can be directly controlled via the magnitude of E. In 
fact the magnitude of even magnetic dipolar interactions 
can be controlled independently of the s-wave interac- 
tions by rotating the magnetic field , but we presume 
that the experimentally easiest option is to use a Fesh- 
bach resonance to adjust a s since only the magnitude of 
the magnetic field needs to be tuned. 

In all the calculations depicted in the following figures 
we limited the minimum value of the scattering length 
to be 17.5ao, where ao is the Bohr radius, corresponding 
to e dd = 0.87. This is because as a s — > the approx- 
imation we made in the calculation of the energy func- 
tional that j3 — 0/Rx be small begins to breakdown (see 
Figures 01 and ■ Furthermore, in the Thomas-Fermi 
approximation used here collapse instabilities can occur 
when £ dd > 1. Indeed, as mentioned previously, be- 
cause of these collapse instabilities one can argue that 
the Thomas- Fermi regime does not exist when e dd > 1, 
see |26|. The maximum value of scattering length was 
set at 600ao, corresponding to e dd = 0.03. 

• Fig. [21 shows the total energy -Etot, as defined by 
Eq. |JBJ, of a condensate with a vortex in an oblate 
trap. In an oblate BEC the dipolar interactions are 
predominantly repulsive, which raises the energy 
relative to the s-wave only case. 

• Fig. |3 plots the radial size R x of a condensate with 
a vortex. In a strongly oblate BEC such as this 
one the majority of atoms interact via the radi- 
ally repulsive part of the dipolar interaction which 
consequently increases R x slightly above the pure 
s-wave value. 
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FIG. 2: Total energy (in the Thomas-Fermi approximation) of 
a dipolar BEC with a vortex in an oblate trap (7 = 5). Solid 
curve: both s-wave and dipolar interactions. Dashed curve: 
s-wave only. The s-wave scattering length a s is measured in 
units of the Bohr radius ao- The energy E to t is measured in 
units of the radial harmonic trapping energy of N atoms. 
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FIG. 3: The radial size R x of a condensate with a vortex in 
an oblate trap (7 = 5). Solid curve: both s-wave and dipolar 
interactions. Dashed curve: s-wave only. R x is measured in 
units of the radial harmonic oscillator length aho of the trap. 
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FIG. 4: The ratio of the vortex core size j3 to the radial size 
R x of a condensate in an oblate trap (7 = 5). Solid curve: 
both s-wave and dipolar interactions. Dashed curve: s-wave 
only. 
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FIG. 5: The aspect ratio k of a condensate with a vortex in 
an oblate trap (7 = 5). Solid curve: both s-wave and dipolar 
interactions. Dashed curve: s-wave only. 



B. Prolate trap 



• Fig. 0] shows (3 = (3/R x , the ratio of the vortex core 
size to the radial size of the condensate. Note that 
this ratio remains small over the chosen range of 
parameters, which, as mentioned above, is neces- 
sary for the self-consistency of the calculation. In 
an oblate trap we see that the effect of the dipolar 
interactions is to increase /3 beyond that found in 
the pure s-wave case. 

• Fig. O depicts the aspect ratio n = R x /R z of a con- 
densate with a vortex. The magnetostriction that 
reduces k by elongating the condensate is clearly 
visible. 

• Fig. El gives the critical angular velocity f2 c of the 
condensate above which it is energetically favorable 
to form a vortex. We see that the effect of dipolar 
interactions in an oblate trap is to decrease Q c . 



We now turn to the case of a prolate trap for which the 
dipolar interactions are predominantly attractive. We 
choose a trap which has the inverse aspect ratio to the 
previous oblate case, namely with frequencies u> x — uj y — 
2ir x 200 rad/s as before, so that the harmonic oscillator 
length of the trap along the x-direction remains the same, 
but with ui z = 2ir x 40 rad/s, so that 7 = 0.2. 

• Fig- D depicts the total energy E tot of the vortex 
state in an prolate trap. The effect of the mainly 
attractive dipolar interactions is to shift the energy 
downwards, oppositely to the case of an oblate trap. 

• Fig. [SI plots the radial size, R x , of a dipolar con- 
densate with a vortex in a prolate trap. R x is 
smaller than in the pure s-wave case because of 
the re-distribution of atoms by the anisotropic in- 
teractions. 

• Fig. EH shows that the ratio of the vortex core size 
to the radial size of the condensate, (3 = (3/R x . 
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FIG. 6: The critical angular velocity £l c above which a vortex 
state is energetically favorable in an oblate trap (7 = 5). 
Solid curve: both s-wave and dipolar interactions. Dashed 
curve: s-wave only. These curves were calculated using Eq. 
J5J. fi c is measured in units of the radial harmonic oscillator 
angular frequency ui x of the trap. 



Note that this ratio remains small over the cho- 
sen range of parameters, which is necessary for the 
self-consistency of the calculation, although not as 
small as in the oblate case. Notice also that in com- 
parison the pure s-wave case f3 is smaller in the 
presence of dipolar interactions in a prolate trap, 
and this trend is opposite to that found in an oblate 
trap. 

Fig, mil gives the aspect ratio k = R x /R z of a con- 
densate with a vortex in a prolate trap. Due to the 
magnetostriction k is smaller than in the pure s- 
wave case. This is the same behavior as was found 
in an oblate trap. 

Fig. II II shows that, in contrast to the oblate case, 
the dipolar interactions in a prolate trap increase 
the value of the critical angular velocity f2 c of the 
condensate above which it is energetically favorable 
to form a vortex. 



IV. AN EXPLICIT EXPRESSION FOR fl c 

The effect that f2 c is decreased in an oblate trap and 
increased in a prolate trap is the main result of this pa- 
per. Let us see if we can interpret this physically. The 
value of f2 c shown in Figures and 1111 for the purely 
s-wave interactions (dashed curve) agrees very well with 
the analytic formula valid in the Thomas- Fermi limit [LH 



2mRl 



h , 0.67R X 
In ■ 



(30) 



This is derived by integrating the kinetic energy den- 
sity n(p, z)v 2 (p)/2, arising from the superfmid flow v s = 
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FIG. 7: Total energy (in the Thomas-Fermi approximation) 
of a dipolar BEC with a vortex in a prolate trap (7 = 0.2). 
Solid curve: both s-wave and dipolar interactions. Dashed 
curve: s-wave only. 
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FIG. 8: The radial size R x of a condensate with a vortex in a 
prolate trap (7 = 0.2). Solid curve: both s-wave and dipolar 
interactions. Dashed curve: s-wave only. 



(H/m)l/p around the vortex, over the profile of the con- 
densate. The lower limit of the integral is set by the 
vortex core size which is given by the healing length £ s 



L = 



1 



\/87m a s 



(31) 



The relevant density appearing in this expression is taken 
to be the density at the centre of the trap no in the 
absence of a vortex, which in the Thomas-Fermi limit is 
given by 



157V 1 
~8^T R Z R 2 ' 



(32) 



In the Thomas-Fermi regime the radii of a condensate 
with s-wave interactions obey 



Ri — a\v 



lhNa s 

flho 



1/5 _ 
£*>h, 



— (33) 



where ZUho = [u x lo v u) z ) x / 3 and — y/h/ (mSJ^J. Com- 
bining all these quantities together one finds that as 
a s — > the critical rotation frequency diverges as f2 c cx 
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FIG. 9: The ratio of the vortex core size fj to the radial size 
R x of a condensate in an prolate trap (7 = 0.2). Solid curve: 
both s-wave and dipolar interactions. Dashed curve: s-wave 
only. 
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FIG. 11: The critical angular velocity £l c above which a vortex 
state is energetically favorable in a prolate trap (7 = 0.2) . 
Solid curve: both s-wave and dipolar interactions. Dashed 
curve: s-wave only. 
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FIG. 10: The aspect ratio k of a condensate with a vortex 
in a prolate trap (7 = 0.2). Solid curve: both s-wave and 
dipolar interactions. Dashed curve: s-wave only. 



— 2/5 

a s loga s , although, of course, eventually the neglected 
quantum pressure term will prevent this divergence. 

Now consider the dipolar case. A formula for fl c closely 
resembling l|3Ufl should still apply. Indeed, the inverted- 
parabola Thomas-Fermi density profile is common to 
both the s-wave and dipolar cases and it is this that 
gives rise to the specific 5/2 and 0.67 numerical factors 
appearing in i|30|) . The one difference we might expect 
concerns the lower limit of the integration which is set 
by the vortex core size. In the presence of dipolar in- 
teractions we can no longer assert that the core size is 
solely determined by £ s since this makes no reference to 
dipolar interactions. Indeed, we might expect that as the 
scattering length vanishes it is replaced by an equivalent 
length scale set by the dipolar interactions of the form 
ad = (Cdd/3)m/(47r?i 2 ). However, expression (I3U|| only 
has logarithmic accuracy |33| and is relatively insensitive 
to the lower cutoff £ s of the kinetic energy integral ap- 
pearing inside the logarithm. The dominant change in 
due to the presence of dipolar interactions is there- 
fore likely to be in the radial size R x . This makes one 
wonder whether the same expression (|3U|1 for £l c approx- 



imately holds for dipolar BECs if the radius R x is modi- 
fied to include the effect of the dipole-dipole interactions 
but the changes in the healin g le ngth are ignored? The 
answer to this question is yes yj] as we shall now show. 

The Thomas-Fermi radii of a vortex-free BEC cylindri- 
cally symmetric in the x — y plane (dipoles aligned al ong 
z) with both contact and dipolar interactions are |25ll26| 

\15gNn f /3k 2 /(k) \1 1 1/5 , . 

and R z = R x /k. The function /(k) appearing in this 
expression is given by Eq. (|18fl and the value of k is de- 
termined by the transcendental Eq. I|22|l. In Fig. I|12|l we 
compare the value of f2 c calculated by the energy mini- 
mization method presented in Section IlIII with the value 
of f2 c obtained from the explicit expression (|30|) where 
R x is given by (|34fl . The agreement is strikingly good, 
with the exception of very small a s in the prolate case. In 
fact, in view of the closeness of the match one is tempted 
to conclude that whilst the long-range dipolar interac- 
tions strongly influence the boundary of the condensate 
and hence the large scales represented by R x , the shorter 
range van der Waals interactions which set the scale for 
a s and hence £ s continue to dominate the shorter range 
physics setting the size of the vortex core, except when 
a s becomes very small indeed. 

V. COMPARISON OF VORTICES WITH 
OTHER TYPES OF ANGULAR MOTION 

We have seen in previous sections that for a dipolar 
BEC in a rotating prolate trap, the formation of a vor- 
tex is increasingly energetically suppressed as the dipo- 
lar interactions grow in strength. Indeed, once f2 c > u) x 
the vortex would be impossible to realize, in a harmonic 
trap at least, because the trap can no longer counter- 
act the centrifugal force and the BEC would fly apart 
(neglecting, of course, the intriguing possibility that the 



9 



„ 0.8 




100 200 300 400 500 600 

a s [ad 



FIG. 12: Comparison of the critical angular velocity Q c cal- 
culated by the methods of Section ITTT1 (solid curves) with the 
prediction of the analytic formula Eq. where the conden- 
sate radius R x is calculated for a dipolar BEC in the Thomas- 
Fermi limit with no vortex (dash-dot curves), see text for de- 
tails. The upper two curves are for the prolate case, 7 = 0.2, 
and the lower two curves are for the oblate case, 7 = 5, as 
before. Note that the match is so good that it is hard to dis- 
cern the difference between the solid and the dash-dot curves, 
except for very small a s in the prolate case. 



attractive nature of the mean-field interactions might be 
capable of holding the condensate together for rotational 
frequencies £1 greater than lu x ). The question then arises, 
how does a prolate dipolar BEC rotate? One possibility 
is for the condensate to generate other types of excitation 
that carry angular momentum, such as quadrupole shape 
oscillations which have an azimuthal angular momentum 
projection m = 2. In the frame of reference rotating 
with the trap a quadrupole oscillation appears as a sta- 
tionary distortion of the density profile in the x — y plane 
|35l l8r1 | and it is generated through a dynamical instabil- 
ity which has a threshold frequency fidyn that in the pure 
s-wave case is typically considerably higher than £l c (in 
the case of an axisymmetric trap Qdyn = i^x/V^). Dy- 
namical instabilities corresponding to shape oscillations 
play an important role in the dynamics of how a vor- 
tex actually enters the condensate, a subject we have not 
touched upon since we have only considered the energet- 
ics rather than the dynamics of vortex formation. In the 
case of pure s-wave condensates it is an experimental fact 
[T^ | that in order to generate vortices the trap must be 
rotated at a frequency considerably exceeding the condi- 
tion of energetic stability given by Eq. I|3U[) and instead 
corresponding closely to fldyn- 

The dynamical instability of a rotating dipolar BEC 
to shape oscillations has recently been investigated [3^ |. 
and it was found that dipolar interactions always lower 
the threshold rotational frequency f^dyn at which the dy- 
namical instability occurs, in both the prolate and oblate 
cases. It is interesting to note that in the pure s-wave 
case the threshold frequency for the instability is inde- 
pendent of interaction strength (even though the insta- 
bility relies on the existence of two-body interactions 0) 
but, conversely, in the dipolar case is strongly affected 



by the magnitude of the interactions (but perhaps, based 
on the discussion of Section II VI the effect of dipolar in- 
teractions is simply through the degree to which these 
change the shape and aspect ratio of the condensate). 
One might tentatively conjecture, therefore, that as the 
critical rotation frequency for energetic favorability of the 
vortex state goes up, and the critical rotation frequency 
for the dynamical instability towards shape oscillations 
goes down, shape oscillations will be formed that will not 
go on to form vortices. 

Let us now also briefly consider another possibility, 
namely oscillations of the centre of mass (CM). In a 
harmonic trap CM oscillations decouple from internal 
excitations provided the interactions are pairwise, and 
therefore take place at the trap frequency. Like low- 
lying shape oscillations, CM oscillations also preserve the 
parabolic density profile and so might be energetically 
preferable to a vortex in the presence of attractive inter- 
actions. Consider a BEC executing a CM oscillation of 
the form 

r cm (i) = r cos(uj x t)x + r a sin(u> x t)y (35) 

where r cm (t) is the position of centre of mass of the BEC, 
which behaves like a classical particle of mass Nm with 
potential energy MUj^ m /2. This circular motion has 
an angular momentum L cm — Nmu> x rQ, with the actual 
value of L cm being determined by the radius rg. Note 
that (|35[) describes dynamics obeying the superfluid irro- 
tationality condition (V x v s = 0) as may be verified by 
observing that the wave function which is an exact so- 
lution of the time-dependent Gross-Pitaevskii equation 
describing this type of motion is of the form [j| 

*(r,i) = exp[-i^t/n]exp{ia(t)[a;- A(t)/2]/ti} 
xcxp{if3(t)[y-B(t)/2]/h} 
x* (x-A(t),y-B(t),z) (36) 

where ^o(x, y, z) is a stationary state satisfying the time- 
independent Gross-Pitaevskii equation, and the parame- 
ters {A(t),a(t);B(t),P(t)} obey a = mA, Amu? x = -a 
and j3 = mB, Brnuj^. — —f3. 

The energy associated with the type of CM oscilla- 
tions given by jSPl is E cm — Nmui x rQ. In order to make 
a comparison between the energy of the vortex state and 
the CM motion we set L cm — Nh, so that the two have 
the same angular momentum. Then r$ = h/(mu! x ), i.e. 
the radius of motion is equal to the trap oscillator length, 
and so E^ = Nhto x . Noting that the ground state 
energy, Eq, is the same for both the vortex and CM mo- 
tion, we can therefore compare E cm with E v = NMl c . In 
order to be competitive with a vortex, CM motion with 
L C m = Nh would therefore require f2 c > lu c , but this is 
exactly the rotational speeds for which the harmonic trap 
no longer confines the atoms and so we surmise that CM 
motion is not energetically favorable over a vortex state 
in a harmonic trap. This does not preclude, however, 
that centre of mass motion takes place in preference to 
other types motion for angular momenta L < Nh. Also, 
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perhaps centre of mass motion can be favored in non- 
harmonic traps. 



VI. CONCLUSIONS 

The question of how quantum fluids rotate is a fas- 
cinating one. The versatility of atomic BECs means 
that this problem can now be studied in systems with 
repulsive or attractive interactions, as well as the case 
of dipolc-dipole interactions which are partially repulsive 
partially attractive and long range. In this paper we have 
used an approach based on an analogy to electrostat- 
ics that allows the explicit calculation of the anisotropic 
mean- field potential inside a dipolar BEC. This approach 
gives quite general insight and in particular shows that 
the long-range part of the interaction is especially im- 
portant in places where the density profile has a large 
curvature in the direction of polarization, such as the 
ends of vortices. Numerical minimization of the total en- 
ergy functional calculated by this method in the Thomas- 
Fermi regime indicates that, in comparison to the case 
of pure s-wave contact interactions, dipolar interactions 
lower the critical rotation frequency f2 c of a BEC neces- 
sary to make a vortex energetically favorable in an oblate 
trap and raise the critical rotation frequency in a prolate 
trap. The same results can also be accurately reproduced 
using the analytic formula Eq. I|30|) well known in the 
usual case of s-wave contact interactions providing the 
modification of the radius R x of the BEC due to dipolar 
interactions is accounted for using Eq. (O. The analytic 
formula allows one to attribute the principal change in 
Q c caused by the dipolar interactions to changes in R x , 
i.e. large scale changes in the overall shape of the con- 
densate, rather than changes on the much smaller scale 
of the healing length which determines the size of the 
vortex core. Finally, we have also discussed angular mo- 
mentum carrying shape oscillations as well as centre of 
mass oscillations as competitors to vortices in rotating 
prolate dipolar BECs, and tentatively conclude that un- 
der certain circumstances shape oscillations such as the 
quadrupole oscillation may be preferred. 
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APPENDIX A: THE ENERGY FUNCTIONALS 

In this appendix we write down the explicit expressions 
for the functionals giving the kinetic, trapping, s-wave 
scattering, and dipolar interaction energy of the conden- 
sate, as defined by Equations Q, ©, (fTU|) . and (fTTf l . 
respectively. Beginning with the kinetic energy, we first 
note that in cylindrical coordinates the gradient is given 
by 



'dp 



ld_ 

' pd<j> 



' dz 



(Al) 




which acts upon the wave function given by the square 
root of the density ansatz (|24l) 



(j),z) 



(A2) 

The phase <f> accounts for the superflow around the vor- 
tex in the usual way, viz. v s How- 
ever, the Thomas-Fermi approximation corresponds to 
neglecting the kinetic energy arising from the curvature 
of the slowly varying condensate background envelope 
— p 2 /R 2 — z 2 /R 2 . We shall consequently drop all 
terms which originate from the gradient of this term. On 
the other hand, it is essential to retain the terms arising 
from the gradient of the vortex part of the wave func- 
tion \J\ — P 2 /(p 2 + (3 2 ) which varies rapidly in the e p 
direction, as well as the superflow term. Under these 
approximations we find the kinetic energy to be 



h n Q irR z 
2m 3 
11/3 4 + 16/3 2 



2G 

y 



+ 



V 7 ! + fi 2 



11/3 2 
— — arctanh 



(A3) 



1/V1 + /3 2 



where no is the number density at the centre of the trap, 
as given by Equation (|2*7I) . 

Using the density ansatz (|24|l the trapping and s-wave 
scattering energies are straightforwardly evaluated to be 
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and 



8tt 2 Rl 
E sw = -gn - 



-p 2 (4 + 13/3 2 + W 4 ) Jl + P 



x arctanh l/yl + (3 



(A5) 



The dipolar energy functional is more difficult to cal- 
culate since the dipolar mean-field potential $dd(r) is 
non-local. In the text we argued that 

E M ^\J d 3 r« bg (r)^(r)+ J d 3 r« v (r)^(r) (A6) 



where $jd( r ) * s §i ven by the quadratic function Equation 
(117(1 . These two integrals can be evaluated explicitly and 
the results are 



\ J d 3 r n bg (r)^(r) = - l -^>e M Rlf(,) (A7) 



and 
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